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We study the deconfining phase transition at nonzero temperature in a SU(N) gauge theory, 
using a matrix model which was analyzed previously at small N. We show that the model is soluble 
at infinite N, and exhibits a Gross- Witten-Wadia transition. In some ways, the deconfining phase 
transition is of first order: at a temperature Td, the Polyakov loop jumps discontinuously from to 
J, and there is a nonzero latent heat ~ iV 2 . In other ways, the transition is of second order: e.g., 

the specific heat diverges as C ~ 1/(T ~ Td) 3 when T — > Tj~ . Other critical exponents satisfy the 
usual scaling relations of a second order phase transition. In the presence of a nonzero background 
field h for the Polyakov loop, there is a phase transition at the temperature Th where the value of 
the loop = 1, with T h < T d . Since dC/dT ~ 1/(T - T h ) 1/2 as T — > T+, this transition is of third 
order. 



The properties of the deconfining phase transition for 
a SU(N) gauge theory at nonzero temperature are of 
fundamental interest. At small AT, this transition can 
only be understood through numerical simulations on the 
lattice [1]. Large N can be studied through numerical 
simulations [2] and in reduced models [3]. In the pure 
glue theory, this transition can be modeled through an 
effective model, such as a matrix model [4-10]. 

One limit in which the theory can be solved analyti- 
cally is by putting it on a sphere of femto-scale dimen- 
sions [11, 12]. An effective theory is constructed directly 
by integrating out all modes with nonzero momentum, 
and gives a matrix model which is soluble at large N [13- 
16]. As a function of temperature, it exhibits a Gross- 
Witten-Wadia transition [17]. That is, it exhibits aspects 
of both first order and second order phase transitions; 
thus it can be termed "critical first order" [12]. Since the 
theory has finite spatial volume, however, there is only 
a true phase transition at infinite N. Thus on a femto- 
sphere, the Gross- Witten-Wadia transition appears to be 
mere curiosity. 

Matrix models have been developed as an effective the- 
ory for deconfinement in four spacetime dimensions (and 
infinite volume). These models, which involve zero [6], 
one [7], and two parameters [8, 9], are soluble analytically 
for two and three colors, and numerically for four or more 
colors. In this paper we show that these models are also 
soluble analytically for infinite N. Most unexpectedly, 
we find that the model exhibits a Gross- Witten-Wadia 
transition, very similar to that on a femtosphere. This is 
surprising because on a femtosphere, the matrix model is 
dominated by the Vandermonde determinant, and looks 
nothing like the matrix models of Rcfs. [6-9]. This sug- 
gests that the Gross- Witten-Wadia transition may not be 
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an artifact of a femtosphere, but might occur for SU (oo) 
gauge theories in infinite volume. At the end of this pa- 
per we estimate how large N must be to see signs of the 
Gross- Witten-Wadia transition at infinite N. 



I. ZERO BACKGROUND FIELD 

We expand about a constant background field for the 
vector potential, 

27rT 

A% = ^j- qi S» , (1) 

where i,j — 1 . . . N. This Aq field is a diagonal SU(N) 
matrix, and so J"'., q^ = 0. The thermal Wilson line is 
the matrix L = exp(27rzq); its trace is the Polyakov loop 
in the fundamental representation, l\ — trL/iV. At any 
JV, this represents a possible ansatz for the region where 
the expectation value of the Polyakov loop is less than 
unity. This region has been termed the "semi" quark 
gluon plasma (QGP) [5]. At infinite N, this ansatz is the 
simplest possible for the master field in the semi-QGP. 
The potential we take is a sum of two terms, 

V cS (q) = - dx(T) Vi(q) + da(T) V 2 (q) , (2) 

where 

K(?)= E \Qi-Qj\ n a-h-Qj\) n - (3) 

The term ~ V 2 {q) is generated perturbatively at one loop 
order; that ~ V±(q) is added to drive the transition to 
the confined phase. Previously, the functions d± and <i 2 
were chosen as d\{T) = (2tt 2 /15) a T 2 Tj and d 2 (T) = 
2tt 2 /3(T 4 - c 2 T 2 Tj) , where T d is the temperature for 
deconfinement [7-9]. These matrix models also included 
terms independent of the q's, ~ C3 T 2 T% and ~ BT^. 
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The values of these parameters were chosen to agree 
with results from numerical simulations on the lattice 
[7-9]. As we show, however, when N is infinite, at the 
transition temperature the nature of the solution is in- 
dependent not only of the values of these parameters, 
but even of the choice of the functions d\ (T) and d 2 (T) 
(modulo modest assumptions, given later). 

The matrix model in Eqs. (2) and (3) is rather different 
from that on a femtosphere [11, 12]. On a femtosphere 
the dominant term driving confinement is the Vander- 
monde determinant, ~ Ilj j log | exp(27rigj) — exp(27ri(?j)|; 
in the present model it is the terms ~ V n {q). The log- 
arithmic singularities of the Vandermonde determinant 
are stronger than those of the from the absolute values 
~ \qi - qj\ n in the ~ V n (q). 

To treat infinite TV, we introduce the variable x = i/N, 
so that qi — > q{x), and the potential is an integral over 
x. It is useful to introduce the eigenvalue density, p{q) = 
dx/dq [13]. The integrals over x then become integrals 
over q, weighted by p(q). The eigenvalue density must 
be positive, and by definition is normalized to 



(4) 



Polyakov loops are traces of powers of the thermal Wilson 
line, 



L = ' tr IP 

3 N 



dq p(q) cos(2tt jq) 



(5) 



As noted before, the first Polyakov loop, £i, is that in the 
fundamental representation. For j > 2, the relationship 
of the tj to Polyakov loops in irreducible representations 
is more involved [5], but all £j are gauge invariant, and 
so physical quantities. 

By a global 0(2) rotation we can assume that the ex- 
pectation value of l\ is real. Consequently, we take p(q) 
to be even in q, p(q) — p(—q). Anticipating the results, 
we also assume that the integral over q does not run the 
full range from — ^ to 1, but only over a limited range, 
from — q to +g - 

Going to integrals over q, we can take out overall fac- 
tors of N 2 from the potentials, with V n (q) = N 2 V n (q), 
where 

V n (q)= [dq f d q 'p(q)p(q')\q-q'\ n (l-\q-q'\r • (6) 



In this expression and henceforth, all integrals over q run 
from —go to +qo, as in Eqs. (4) and (5). 

We then define V e g(q) = N 2 V c s{q), where Kff = 
— dxVi + d 2 V 2 . Solving the model at infinite N, then, 
is just a matter of finding the (minimal) stationary point 
of V e f[(q) with respect to the q^s. 

The equations of motion follow by differentiating the 
potential in Eq. (2) with respect to qt, and then taking 
the large N limit. Doing so, we find 



= [di + d 2 ] q - 
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+ d 2 j dq'p(q') [-3(q - q')\q - q'\ + 2(q - q'f] , (7) 

where sign(a;) = ±1 for x ^ 0. For simplicity we write 
d\(T) and d 2 (T) just as d\ and d 2 . 

To solve the equation of motion in Eq. (7), we follow 
Jurkiewicz and Zalewski [16] and use the following trick. 
What is difficult is that Eq. (7) is an integral equation 
for p(q). To reduce this to a differential equation, take 
<9/<9gofEq. (7), 



= d 1 



di p(q) 



(8) 



+6d 2 J dq / p(q / )[-(q-q / )sign(q-q r ) + (q-q / ) 2 ] . 

Notice that this does not give us the second variation of 
the potential with respect to an arbitrary variation of q, 
which is related to the mass squared. Instead, we take 
the derivative of the equation of motion, with respect to 
a solution of the same. 

We then continue until we eliminate any integral over 
q' . Taking d/dq of Eq. (8) gives 



di^=6d 2 j dq'p { ,/) -sv'uUi .,/.)-:■ -1U, 



Lastly, by taking one final derivative, we obtain 
d 2 



2 p(q) + d 2 [p(q)-l} =0 



dq 



In this expression we introduce the ratio 

AT) - ^ . 
di(T) 



(9) 
(10) 

(11) 



dq' p(q') sign(q - q 1 ) 



We assume that like the solution at small N [7-9], that 
d(T) increases with T, and d(T) — > oo as T — >• oo. We 
note that the only detailed property of d(T) which we 
require is that its expansion about Td is linear in T — Tj. 
This is a minimal assumption which is standard in mean 
field theory. 

We thus need to solve Eqs. (7) - (10), subject to the 
condition of Eq. (4). The solution of Eq. (10) is trivial, 

P(q) = 1 + bcos(dq) , q:-q ^q Q , (12) 

where b is a constant to be determined. We assume that 
p(q) = for \q\ > qg. We have checked numerically that a 
multi gap solution [16], where p(q) ^ over a set of gaps 
in q, does not minimize the potential; see the discussion 
at the end of Sec. (III). 

When q < |, p(qo) 7^ 0, and the solution drops discon- 
tinuously to zero at the endpoints. This stepwise discon- 
tinuity is charactertistic of the model, and presumably 
reflects the singularities from the absolute values in the 
potential. 
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The eigenvalue density in Eq. (12) is simpler than that 
in the Gross- Witten model [11, 12, 14-16], where 



Pow(q) = ^cos(Trq) 



sin 2 (7rg) 
sin 2 (7rg ) 



1/2 



(13) 



For any g , this vanishes at the endpoints, Pawi^Qo) = 
0, while at the transition, go = \. Due to the Vander- 
monde determinant in the potential, the density pgw(q) 
has a nontrivial analytic structure in the complex q- 
plane, while p(q) does not. Since the Vandermonde po- 
tential is so different from V c g, though, it is natural to 
find that pgw(q) is unlike p(q) in Eq. (12). 

Eq. (12) solves Eq. (9) without further constraint. To 
solve the remaining equations, remember that all inte- 
grals run from —go qo- The normalization condition 
of Eq. (4) gives bsm(dq ) = d(| — qo)- After some alge- 
bra, one can show that Eqs. (7) and (8) are equivalent, 
with the solution 



cot (d qo) = 



d [I 



and 



b 2 = 



3 V 2 



'7o 



'7o 



1 



d(l/2-®,) 
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'7o 



(14) 



(15) 



Thus in the end, we only have to solve two coupled al- 
gebraic equations, Eqs. (14) and (15), for g and b as 
functions of d = d(T). 

At low temperature, d is small, and the theory is in 
the confined phase, where b = and go = \ ■ The eigen- 
value density is constant, p(q) = 1, and all Polyakov 
loops vanish, lj = 0. Thus the confined phase is char- 
acterized by the maximal repulsion of eigenvalues. The 
Gross- Witten model also has a constant eigenvalue den- 
sity in the confined phase, which is expected, as only a 
constant eigenvalue density gives lj = for all loops. 

In the limit of high temperature d — ¥ oo. The solution 
is go — 6/d 2 and b = d 2 /12. The eigenvalue density 
is p w gP/12, which becomes a delta-function S(q) for 
infinite d. That is, at high temperatures all eigenvalues 
coalesce into the origin, and all Polyakov loops equal one, 

l i = l - 

As the temperature and so d(T) is lowered, the tran- 
sition occurs when g = |, for which d(T^) = 2ir. At the 
transition point, the eigenvalue density is 

p(q) = 1 + cos(27rg) ; T = T d . (16) 

From Eq. (5), 

Zi(T+) = \ , ^(T d ) = 0,j>2. (17) 

Thus at the transition, only the Polyakov loop in the 
fundamental representation is nonzero, equal to \. 

What is unforeseen is that at T£ , the eigenvalue den- 
sity in the present model, Eq. (16), coincides identically 



with that in the Gross- Witten model, Eq. (13). Conse- 
quently, properties exactly at T~[~ , such as the expectation 
values of the lj, are the same in the two models. Since 
they differ away from T d , other properties are similar, 
but not necessarily identical. 

Consider the behavior in the deconfined phase just 
above the transition point, taking d — 2tt(1 + Sd). The 
solution is g^ = 5(1 — Sq), where 



(18) 



6=1 



1 /25tt 2 



2\ 1/5 



Sd 2/5 - 



29 /25tt 2 



2\ 2 /5 



2 V 3 J 56 V 3 

Using this, one finds that 



Sd i/5 



+ . 



1 1 /25tt 



2 4 V 3 



2 \ 1/5 



Sd 2/5 



(19) 



(20) 



while all lj ~ Sd for j > 2. 

Remember that at T^, i\ jumps discontinuously, from 
to I , as expected for a first order transition. Assuming 
that 5d~T d - T, though, Eq. (20) shows that as T -> 
T+ 



h(T) 



(T d - Tf 



(21) 



That is, near the transition £i(T) exhibits a power like 
behavior which is characteristic of a second order phase 
transition — although £i(T£) ^ 0. 

For arbitrary d, after some algebra one finds that at g^ , 
the solution of Eqs. (14) and (15), the potential equals 



vm)-v^ 



conf 



15 V 2 



(22) 



The potential in the confined phase is V^g nf = V c s(^) = 
—di/6 + <i2/30. In these matrix models, the pressure is 



P(T) = -V BS (q S o) + Kff 



(23) 



This subtraction ensures that the pressure, and the asso- 
ciated energy density, are suppressed by ~ 1/7V 2 in the 
confined phase. In the models of Ref. [7, 8], V^g is 
given by the term ~ C3. Expanding about T d , 



V eS (qo)-V c ^ 



48d 2 



Sd 



270d 2 
7tt 3 



95 \ 1/5 

5) M " B - 



(24) 

Assuming that Sd ~ T — T d , as is true of the functions 
in Refs. [7-9], the leading term in Eq. (24) ~ Sd shows 
that the first derivative of the pressure with respect to 
temperature, which is related to the energy density e(T), 
is nonzero at T d . Since the pressure and the energy 
density are suppressed by ~ 1/N 2 in the confined phase, 
the latent heat is nonzero and ~ TV 2 , ~ e(Tf). 
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Using the explicit forms for d\(T) and d,2(T), we find 
that the latent heat is e(T+)/(7V 2 T|) = 1/tt 2 - .10... 
This is about four times smaller than the lattice results 
of Ref. [2] who find ~ 0.39 for the same quantity. The 
lattice results can be accomodated by adding a term like 
a MIT bag constant to the model [8]. Such a term is 
~ Tj but independent of the q's, and so only changes the 
latent heat, but does not affect any other result. 

The second term in Eq. (24) shows that the second 
derivative of the pressure with respect to temperature 
diverges as T - 



d 2 
dT 2 



P(T) 



1 



(T - T d )« 



(25) 



This is the usual divergence of the specific heat for a 
second order phase transition. 



II. NONZERO BACKGROUND FIELD, T = T d 



and 

2 V 12 J 56 V 80 



4\ V5 



(30) 

For this solution, at the minimum the /i-dependence of 
the potential is 



di 



di 



^0^) = -^+ 112 , 127r 3 



25 



1/5 



The expectation value of the loop l\ is 
2 4 V 12 / 112 V 80 



(31) 



(32) 



Hence l\ — \ ~ h 1 / 8 , where 5 — 5/2. This shows that the 
critical exponents of this model satisfy the usual Griffths 
scaling relation, 



Background fields can be added for each loop £j. In 
this paper we just consider a background field for the 
simplest loop, t\, since only that is nonzero at Td, Eq. 
(17). We add 



V h {q) 



di 



h£i 



(26) 



to the potential V g(q), and find the solution as before. 
After taking three derivatives of the equation of motion, 
with respect to a solution, we obtain the analogy of Eq. 
(10), 



dq- 



: p(q) + d 2 [p(q) - 1] + (2tt) 2 h COs(27rq) = . (27) 



This equation is valid for any d. It is necessary to treat 
the case of T d , where d = 2n, seperately from T ^T d . 

In this section we consider the point of phase transi- 
tion, where d = 2ir. The solution of Eq. (27) is 



p(q) = 1 + 6 cos(27r q) — nh q sin(27rg) 



(28) 



where q : —q^ — > qo. Notice that the /i-dependent term 
q sin(27T(7) arises because when T = Td, Eq. (27) repre- 
sents a driven oscillator at the resonance frequency. The 
value of the constants 6 and go now depend upon both 
d(T) and the background field, h. 

The analogy of Eq. (9) is solved by Eq. (28). The 
normalization condition, Eq. (4), plus the analogy of Eq. 
(8), gives two equations for 6 and qo; as before, Eq. (7) 
does not give a new condition. 

When ft, ^ 0, the explicit form of the analogy of Eq. 
(4) is elementary, but that of Eq. (8) is rather ungainly. 
We thus present the results of the solution in the limit 
of small background field, h <C 1. We find that q^ — 



(1 — Sq), where 



Sq = 



' 45 
2tt 4 



1/5 



3/3 
14 I 200 ir 2 



1/5 



(29) 



2 - a = 0(1 + S) . 



(33) 



The effective potential, as a function of £% , is computed 
by taking the Legendre transform, 



( 27r ) 

Expanding the potential in 5t\ = l\ — % at , 



128^3(^1 

35 7T 3 



6£[ 



7/2 



32 d x 



(34) 



(35) 



This is a very flat potential, starting only as (l\ — \) 7 ^ 2 - 
This is in contrast to the femtosphere, where the po- 
tential behaves as ~ {l\ — |) 3 about the similar point 
[11,12]. 

Expanding at T d gives the expansion of the potential 
about £\ = 0. One can show, and we verify in the next 
section, that this potential vanishes. This implies that 
the potential has an unusual form: it is zero from l\ : 
— > 1, and then turns on as in Eq. (35). Graphically, 
this potential is like that on the femtosphere; see, e.g., 
Fig. (1) of Ref. [12]. 



III. NONZERO BACKGROUND FIELD, T T d 

Consider now the theory in a nonzero background field 
for £i, Eq. (26), away from the transition, so d ^ 2ir. The 
eigenvalue density again solves Eq. (27). The solution is 
simpler when d ^ 2-7T, and is just the sum of the solution 
for h — and an ^.-dependent term, 



p(q) = 1 + 6 cos(rfg) + 



1 



1 - (d/2ny 



h cos(27rg) . (36) 



The solution follows as previously, and we simply sum- 
marize the results. 
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We first consider the confined phase, defined to be the 
solution for which qo — \ and 6 = 0. The expectation 
value of the loop l\ is 



1 h 



1 l-(d/2n) 2 2 
For this solution the potential equals 



(37) 



h 2 



Performing the Legendre transformation, we find 



(39) 



This shows that in the confined phase, when d < 2ir the 
mass squared of the i\ loop is positive, as expected. It 
also shows that this mass vanishes at Td when h = 0; this 
justifies the statements about the potential at the end of 
the previous section. 

Consider a special value of d, d\ = 4ir 2 (l — h); the 
corresponding temperature is defined to be Th, d(Th) — 
dh- At this temperature, the eigenvalue density of Eq. 
(36) coincides exactly with that at the transition in zero 
background field, Eq. (16). Notably, the values of the 
loop at h and T = Th are the same as for h = and 
T = T d : h(T h ) = I with = for j > 2, Eq. (17). 
Thus we may suspect that something special happens at 
T = Th. For example, the confined phase is only an 
acceptable solution when T < Th, as only then is the 
eigenvalue density positive definite. 

This suggests that a phase transition occurs at dh . To 
show this, we compute for about this value of d, taking 
d 2 = d\ + 4-7T 2 h 5d. Solving the model as before, in the 
deconfined phase the solution is — |(1 — ^l): where 



Sq = 



1/2 



W 2 



V6 
40tt 



(8/i - 5) <M 3/2 



b = — (1 - h) 3/2 csc(VT^h ?r) 5d 5/2 



(40) 



(41) 



With this results we compute the potential in the decon- 
fined phase, to find 



v cS (h)-v:rw 



3^6 

'5^3 



8d 5 > 2 



(42) 



Taking Sd ~ Th — T, we find that the third derivative 
of the pressure, with respect to temperature, diverges at 



<9 3 1 

df* P(T) ~ (T — Th) 1 / 2 1 T ^ T t 



(43) 



In zero background field, then, there is a critical first 
order transition at a temperature Td- Turning on a back- 
ground field ~ ht\, the first order transition is immedi- 
ately wiped out for any h =/= 0. Even so, there remains a 



third order phase transition, at a temperature Th < Trf, 
where the expectation value of the loop l\ = h. This 
behavior is the same as on a femtosphcre [11, 12]. 

In principle one can also add a background field for 
any loop, £j for j > 2. It is direct to derive the equa- 
tions of motion and obtain a solution for the eigenvalue 
density. Obtaining the minimum of the potential is not 
elementary, though. The original model of Gross and 
Witten [14] involves the Vandermonde determinant plus 
a term ~ |trL| 2 . The solution for the eigenvalue den- 
sity is a function which is nonzero on one interval, be- 
tween —qo and go- Jurkiewicz and Zalewski [16] showed 
that when terms such as |trL 2 | 2 are added to the Gross- 
Witten model, that in general it involves functions which 
are nonzero on more than one interval. We have checked 
numerically that when only h% ^ 0, that such multi-gap 
solutions do not minimize the potential. We do find, how- 
ever, that multi-gap solutions do minimize the potential 
in the presence of background fields for £j when j > 2. 
Since only l\ ^ at Td and Th, we defer the problem of 
background £j for j > 2. 



IV. FINITE N 

The model can be solved numerically at finite N. This 
confirms, as expected on general grounds [8], that the 
deconfining transition is of first order for any N > 3. It 
also shows that the critical behavior found at infinite N 
is smoothed out at large but finite N. 

Using the numerical solution of the model, in the Fig- 
ure we show the behavior of the specific heat, divided by 
N 2 — 1, for different values of N. To see the putative 
divergence of the specific heat at infinite N, rather large 
values of N are necessary, N > 40. 

This Figure also shows that the increase in the specific 
heat only manifests itself very close to the transition, 
within ~ 0.2% of T d - At present, direct numerical simu- 
lations on the lattice treat moderate values of N ~ 4— 10 
[2]. For most quantities there seems to be a weak varia- 
tion with TV. 

The present matrix model suggests that very near T d , 
a novel phase transition may occur at large TV. The val- 
ues of iV at which critical first order behavior arise can 
presumably be studied only in reduced models [3]. 

This begs the question of whether or not the Gross- 
Witten-Wadia transition does in fact occur at infinite N 
[18]. On the femtosphere, one can easily solve the model 
in the presence of additional couplings, such as (|trL| 2 ) 2 . 
Such additional couplings turn the Gross- Witten- Wadia 
transition into an ordinary first order transition [12]. We 
have not been able to solve the present model in the 
presence of additional couplings. 

The most likely possibility is that as on a femto- 
sphere, the presence of additional couplings washes out 
the Gross- Witten- Wadia transition. Nevertheless, gauge 
theories are remarkable things. Certainly it is worth 
studying SU(N) gauge theories, at very large values of 



() 




N, to see if there is a Gross- Witten-Wadia transition at 
infinite N. 



1.002 



FIG. 1. Plot of the specific heat, divided by (N 2 - 1)T 3 , for 
different values of N. 
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